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I.  INTRODUCTION 


\ 

The  purpose  of  the  study  presented  herein  is  to  demonstrate  the 
applicability  of  a  phenomenological  relation  between  stress,  strain  rate  and 
temperature  in  describing  the  behavior  of  viscoelastic  materials,  in  general, 
and  solid  propellants  in  particular.  The  basic  relation  and  its  specific 
equations  associated  with  various  testing  conditions  are  discussed  in  the 
appendix.  ( 

The  procedure  was  to  extract  from  the  literature  those  test  data 
which  would  not  only  provide  the  required  material  constants,  but  at  the 
same  time  permit  an  evaluation  of  the  relation's  range  of  applicability. 

Although  the  majority  of  the  data  uncovered  left  much  to  be  desired 
in  the  way  of  completeness  (from  the  viewpoint  of  utilization  in  this  evaluation), 
nevertheless.it  was  possible  to  detect  certain  encouraging  trends.  The 
manner  in  which  the  available  data  wereutilized  and  the  interpretive  con¬ 


clusions  are  presented. 
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II.  STRESS-STRAIN  (CONSTANT  STRAIN  RATE) 

The  majority  of  published  experimental  information  for  visco¬ 
elastic  materials  is  in  the  form  of  stress- strain,  ultimate  stress,  or  ultimate 

strain,  constant  strain  rate  data.  In  order  to  extract  from  these  types  of 

AH 

data  the  governing  material  parameters  <r0.  -g-»-  and  s  (or  loge  s)  it  is 
necessary  to  initially  assume  that  Equation  (A5),  with  its  inherent  limitations 
regarding  the  values  of  (J  and  M,  applies  when  the  strains  are  "large.  "  Two 
methods  of  extracting  the  constants  from  the  data  were  derived. 

The  first  method  utilizes  Equation  (A5)  in  the  form 

~  +  lo8e  8  *  “  =  loSe  Y  U> 

„  .....  ...  1  .  AH 

For  any  one  test  condition,  this  equation  is  linear  in  — ,  loge  s  and  -=— , 

(Tq  K 

with  constant  coefficients.  For  any  three  test  conditions  of  a  material 

yielding  stress  data  at  high  strain,  we  have  a  system  of  three  linear  equations 

in  three  unknowns.  The  determinant  of  the  system  is 
I  1  ! 


1 

'  T  l 

1 

1 

‘  t2 

1 

1 

crf  3 

' 

A  =  13  '  Jl  +  °*1  '  *3  ,  °I2  '  °fl 
Ti  T2  T3 


(2) 
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Data  at  no  leas  than  two  temperatures  are  necessary,  since  when 
Tl  =  T2  =  T3,  A  =  0,  (i.  e.  ,  the  system  is  inconsistent  and  has  no  solution) 

~  ,  1  .  j  AH 

Solving  for  — ,  log.,  s  and  — ,  we  have 

On  C  K 
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This  approach  has  the  advantage  of  computing  the  material  constants 
independently  and  utilizing  a  minimum  of  the  test  conditions  and  data  in  a 
straightforward  procedure.  However,  due  to  the  manner  in  which  much  of 
the  data  was  presented  in  the  literature,  a  second  method  was  found  to  be 
more  advantageous  for  this  study. 

If,  for  any  one  material,  there  exists  stress  data  at  large  strains 
acquired  from  tests  at  two  or  more  constant  strain  rates  (€j,  e2,  .  .  .  e^) 
at  the  same  constant  temperature  (Tj  =  T2  =  .  .  .  Tj),  then  Equation  (A5) 
provides  i  -  1  relations  for  the  material  constant  <r0  of  the  following  form 


ffli  ~  **1 

lo8e  T2- 

«1 


(6) 


where  tru  is  the  lowest  limiting  stress  in  the  group. 

AH 

Knowing  the  value  of  rQ  ,  it  then  becomes  possible  to  solve  for  — , 

R 

provided  values  of  the  limiting  stress  exist  for  one  strain  rate  =  e2  =  ...  e.)j 


at  two  or  more  temperature*  (Tj,  T .  .  .  Tj).  Again,  from  Equation (A5), 
we  have  i  -  1  relations  of  the  type 
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9H  -  '<1  ,  ,  T1 

-  +  lo8e  r- 

o-o  Tt 

_1_  _1_ 

Ti  ‘  Ti 


Finally,  the  value  of  s  may  be  acquired  from  the  relation 
e: 


a  = 


Ti  exp 


°ji  _  AH 

o-o  ' RT; 


(7) 


(8) 


where  <rQ  and  are  the  values  acquired  from  Equations  (6)  and  (7), 
respectively  and  e j,  Tj,  and  ir^  are  the  test  data  from  anyone  of  the  previously 
used  constant  strain-rate  tests. 

Theoretically,  in  Equations  (6)  and  (7),  only  two  values  of  <rj  are 

required;  in  Equation  (8),  only  one  value  is  necessary.  However,  if  more 

values  exist,  these  may  be  used  to  provide  an  arithmetic  average  in  the 

AH 

in  the  event  that  minor  deviations  for  <r0,  -j^-,  *nd  8  arise  as  a  result  of 

slight  irregularities  in  the  recorded  or  plotted  data  and  to  also  verify  the 
fact  that  the  parameter  is  indeed  constant  over  a  range  of  strain  rates. 

Mechanical  properties  of  various  propellants  are  presented  by 
Burgardt,  Frazer  and  Britton^ )  for  temperatures  ranging  from  233"  to 
311*K,  and  loading  rates  of  0.  05,  0.2,  2.0,  lO.Oin/min.  A  strain  gage 
length  of  2  inches  was  assumed  to  obtain  strain  rates.  Two  propellants 


were  selected  as  examples. 
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ah 

In  deriving  the  parameter s  ff0,  and  loge  s  for  Propellant  1,  the 
ultimate  stress  data  for  all  loading  rates  at  277.  4°K  were  used,  along  with 
the  stress  at  298®K,  loading  rate  2.  0  in/min. 

The  ultimate  stresses  at  other  temperatures  were  then  predicted 
by  again  using  Equation  (1),  i.  e. 


AH 

RT 


+  loge  e  -  loge  T  -  loge  s 


where  <rD,  and  loge  s  were  the  previously  determined  constants  and  <rj 
was  calculated  for  each  T  and  e.  The  computed  values  compare  favorably 
with  the  experimental  values  given  for  other  temperatures.  The  ratios  of 
experimential  to  calculated  ultimate  stresses  are  given  in  Table  1. 

The  same  method  was  employed  to  compute  the  material  constants 
of  the  second  propellant;  however,  the  ultimate  stresses  used  were  obtained 
at  298°Kfor  all  loading  rates  and  at  277. 4'F,  loading  rate  2.0  in/min. 
Ultimate  stresses  were  then  calculated  and  Table  2  lists  the  ratios  of  experi¬ 
mental  to  calculated  stresses  for  the  second  propellant. 

The  calculated  stresses  are  consistent  with  the  experimental  data 
for  both  propellants  with  the  possible  exception  of  data  at  233'  K.  A  reason¬ 
able  estimate  of  the  limiting  temperature  of  Equation  (1),  is  one-half  the 
absolute  melting  temperature. 

In  the  study  of  Smith's  data ^  for  polyisobutylene  polymers  the 
process  was  again  utilized  and  the  material  constants  were  found  to  be 
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TABLE  1.  RATIO  OF  EXPERIMENTAL  TO  CALCULATED 
ULTIMATE  STRESSES  (PROPELLANT  1) 


Temperature 

Rate  of  Loading, 

°K 

°F 

in/ min 

Ratio 

310.  8 

100 

0.  05 

1.  173 

0.  2 

.  889 

2 

1.  085 

10 

.  923 

298 

77 

0.  05 

1.  173 

0.  2 

1.  188 

2 

1.  000 

10 

1  164 

277.  4 

40 

0.  05 

.  977 

0.  2 

.  980 

2 

•  989 

10 

1.  043 

266.  3 

20 

0.  05 

1.  098 

0.  2 

1.  135 

2 

1.  073 

10 

1.  074 

255.  2 

0 

0.  05 

.  985 

0.  2 

1.  119 

2 

1  008 

10 

1.  051 

233. 

-40 

0.  5 

1.  196 

0.  2 

1.  179 

2 

1.  178 

10 

1.  441 
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TABLE  2.  RATIO  OF  EXPERIMENTAL  TO  CALCULATED 

ULTIMATE  STRESSES  (PROPELLANT  2) 

Temj 

aerature 

Rate  of  Loading, 

•K 

*  F 

in/ min 

Ratio 

310.  8 

100 

0.  05 

1.  046 

0.  2 

.  920 

2.  0 

1.044 

10.  0 

-  921 

298 

77 

0.  05 

1.  027 

0.  2 

1.  000 

2.  0 

1. 010 

10.  0 

1.  101 

277.  4 

40 

0.  05 

.  846 

0.  2 

.  930 

2.  0 

1.  000 

10.  0 

1.  040 

266.  3 

20 

0.  05 

.  970 

0.  2 

1.  046 

2.  0 

1.  160 

10.  0 

1.  185 

255.  2 

0 

0.  05 

1. 081 

0.  2 

1.  083 

2.  0 

1.  222 

10.  0 

1.  379 

233 

-40 

0.  05 

1.  259 

0.  2 

1.  000 

2.  0 

1.  478 

10.  0 

1.  716 
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<rQ  =  3.7,  =  5200,  s  =  0.  37.  The  ultimate  stresses  were  then  calculated 

and  can  be  found  in  Table  3  and  Figure  1.  The  experimental  values  uaed  to 
compute  the  parameters  are  noted  in  the  table. 

In  addition  to  the  ultimate  strain  data,  experimental  stress  -  strain 
data  were  also  given.  Therefore,  the  constants  were  utilized  in  Equation 
(A2)  and  stress- strain  curves  were  generated  for  four  constant  strain  rates 
at  T  =  358* K.  These  curves  appear  in  Figure  2.  The  difference  in  the  high 
strain  rate  stress- strain  curves  is  due  to  the  initial  hypotheses  in  which  the 
steady- state  creep  condition  was  considered. 

The  material  constants  for  a  GR  S  rubber  were  calculated  from 
data  acquired  by  Smith^)  at  temperatures  298*,  316°,  366° K  and  found  to 

Ah 

be  <r0  =  15,  -g—  =  9430,  and  s  =  5  5  60.  The  ultimate  stresses  for  several 
strain  rates  at  360.  8*K  were  computed  and  compared  with  Smith's  experi¬ 
mental  data  in  Figure  3.  Smith's  data  are  based  on  averages  of  two  tests 


at  each  strain  rate. 


TABLE  3.  ULTIMATE  STRESS  OF  POLYISOBUTYLENE 


Strain 

Experimental 

Calcula 

Temperature 

Rate, 

Ultimate  Stress. 

Ultimate 

*K 

•c 

in/in/ sec 

psi 

psi 

241. 4 

-31.  6 

. 000104 

34.  6 

31. 7 

. 00104 

41.  6 

40  9 

.0104 

46.  8 

50.  1 

.  104 

74.  1 

59.  3 

111 

-  1.0 

. 000104 

_  _  _  _ 

21. 4 

. 00104 

29.  3 

30.  6 

.0104 

39.  5 

39-  8 

.  104 

48.  6 

49.  0 

300 

26.  7 

. 000104 

12.  8* 

13.  8 

. 00104 

24.  3* 

23.  0 

.0104 

32.  6 

32.  2 

.  104 

41.4 

41. 4 

333 

60 

.000104 

5.  Of 

6.  6 

. 00104 

14.  3* 

15.  8 

.  0104 

25.  1* 

25.  0 

.  104 

30.  4 

34.  2 

358 

85 

. 000104 

2.  8f 

1.9 

. 00104 

9.  9* 

11  1 

.  0104 

17.  7f 

20.  3 

.  104 

27.  9* 

29.  5 

*Uaed  in  calculating  parameters. 


FIGURE  1 .  ULTIMATE  STRESS  OF  POLYISOBUTYLENE  AT 
VARIOUS  TEMPERATURES  AND  STRAIN  RATES 


FIGURE  2.  STRESS -STRAIN  CURVES  FOR 
POLYISOBUTYLENE  AT  358°  K 


13 


!  »*/  Frm  r/u.M'  trr 


FIGURE  3.  ULTIMATE  STRESS  OF  A  GR-S  RUBBER  AT  360.  8’K 
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III.  STRESS  RELAXATION 


Although  data  are  limited  for  the  teat  conditions  where  strain  and 

temperature  are  constant  and  the  stress  varies  with  time,  Tobolsky^ 

has  published  applicable  data  for  unfractionated  polyisobutylene  at  50% 

elongation  at  temperatures  303®,  323®,  343°,  373°K. 

AH 

To  obtain  the  constants  vQ,  -g— ,  and  s  we  first  utilize  Equation  (A15) 


tanh  •£—  ~  tanh  eXp  (  .  -L  \ 

2*o  2<r0  l>\  t<  ) 

Since  t'  is  dependent  upon  the  temperature  and  contains  two  unknown  param - 
eters,  we  must  separately  obtain  tk  and  t{>.  It  was  necessary,  therefore,  to 
choose  at  T  -  Ta  -  373°K  two  stresses  and  their  times  of  occurrence  which 
gave  a  system  of  two  equations,  i.  e. , 


*1 

tanh  - —  =  tanh 
2<r0 


(9) 


*2  *(t0>  /  t2  \ 

tanh -  -  tanh  — -  exp  l  — —  I 

2*o  2*o  '  la  / 

where  <r  (tQ)  is  the  known  initial  stress  and  tj  -  IX  10*2,  t^  1  1  X  10**. 
Solution  of  »his  system  provides  cr0  and  ta.  Solution  of  Equation  (9)  was  then 

I 

repeated  at  T  *  Tb  -  323° K  to  obtain  tb-  The  applied  stress  <r(t0)  must  be 
known  for  each  test  condition.  In  this  study,  <r(tQ)  was  assumed  to  be  very 
near  the  stress  at  t  =  1  X  10"3;  this  was  proven  to  be  a  good  approximation. 


The  known  value*  oft'  are  used  to  generate  from  Equation  (A16) 


a  new  system  of  simultaneous  equations 


<r0  exp 


*a  = 


*b 


'  RT  a ' 


2s  Ta  Ea 
(Tq  exp 


V  RTb  / 


(10) 


2s  Tb  Eb 


where  E  =  - ~ ~  and  is  dependent  on  the  respective  temperatures.  We 

ah  5 

found  or 0  ■  4,  — —  =  7040,  s  =  2.93  X  10  .  Stress  decay  curves  were  then 
calculated  by  use  of  Equation  (A15)  and  are  shown  in  Figure  4.  The  exper 
mental  and  calculated  curves  are  in  good  agreement. 


&  7-/3.M&S  -P*' 


1 


FIGURE  4.  STRESS  DECAY  OF  UNFRACTIONATED 
POLYISOBUTYLENE  AT  50%  ELONGATION 
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IV.  DISCUSSION 

The  degree  of  verification  of  the  basic,  phenomenological  relation 

for  viscoelastic  materials  is  necessarily  restricted  by  the  rather  limited 

amount  of  applicable  data.  A  brief  literature  study  was  conducted  but  much 

of  the  data  uncovered  was,  for  our  purposes,  incomplete.  From  the  data 

found  and  presented  herein,  however,  certain  encouraging  trends  are  seen. 

AH 

Procedures  for  deriving  the  materials  constants  <r0,  — — ,  and  s 

K 

are  adaptable  to  the  amount  of  usable  data  available.  With  the  restriction 
to  temperatures  at  which  the  viscoelastic  properties  are  in  effect,  ultimate 
stresses  can  be  predicted  within  a  close  percentage.  The  stress-strain 
curves  calculated  correlate  readily  with  the  experimental  curves,  especially 
at  slow  strain  rates.  The  calculated  stress  decay  curves  agree  both  in 

magnitude  and  in  general  form  with  the  experimental  curve,  even  though  the 

\ 

actual  initial  load  was  not  known  and  had  to  be  estimated. 

Within  the  limits  of  this  study,  these  results  are  favorable,  and 
the  application  of  the  basic  phenomenological  relation  to  other  testing 
conditions  for  viscoelastic  materials  should  be  no  less  encouraging. 
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APPENDIX 

The  following  relations  were  obtained  from  Stowell'sl  ^1)  develop¬ 
ment  and  are  repeated  here  for  convenience. 

€  =  Strain 

e  =  Strain  rate 

t  =  Time 

<r  =  Stress 

T  =  Temperature,  0 K  (unless  otherwise  indicated) 

£  =  Young's  modulus 

a  =  Linear  expansion  coefficient 

AH  s  Activation  energy,  cal  per  mole 

R  =  Gas  constart,  taken  as  2  cal  per  mole  per  "K 

8  =  constant 

o'o  =  constant 
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PHENOMENOLOGICAL  CONCEPTS  FOR 
VISCOELASTIC  MATERIALS 


A.  General 

The  basic  phenomenological  relation,  which  forms  the  foundation 
for  mathematically  describing  the  behavior  of  solid  propellants,  received 
its  initial  impetus  from  a  paper  by  Kauzmann(Al).  in  this  paper,  Kauzmann 
following  a  suggestion  of  Condon,  attempted  to  compute  the  creep  rate  for 
a  metal  under  steady  stress  from  the  viewpoint  of  reaction- rate  theory  as 
used  by  physical  chemists,  notably  Eyring.  Although  the  Kauzmann  equation 
contained  much  truth,  experimental  evidence  clearly  indicated  modifica¬ 
tions  were  necessary.  Three  modifications  plus  the  inclusion  of  terms  for 
the  effects  of  elasticity  (always  present  to  some  extent)  and  thermal  expan¬ 
sion  led  to  the  basic  formulation  by  Stowell^^.  This  relation  has  since 
been  used  successfully  to  predict  the  behavior  of  a  variety  of  viscoelastic 
material  under  various  time  and  temperature  conditions  (for  example, 
see  Ref.  A3  and  A6). 

The  basic  equation  proposed  for  a  viscoelastic  material  is  given  by: 


de 

dt 


(Al) 


Essentially,  the  preceding  phenomenological  relation  requires  only  the 
acquisition  of  the  material's  constants  (AH,  s,  and  <rQ)  from  any  one  of 
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the  several,  more  familiar  quasi- static  tests.  Solutions  to  problems 
involving  viscoelastic  materials  are'  obtained  as  integrals  of  Equation  (Al) 
subject  to  the  appropriate  load,  deformation  and  temperature  conditions 
(see  below)  Thus,  there  is  no  recognizable  "equation  of  state,"  although 
Equation  (Al)  may  be  appropriately  called  an  "equation  of  rate  of  change  of 
state.  " 


B.  Stress -Strain  Relations 

Equation  (Al)  yields  the  following  stress-strain  relation  for  constant 
temperatures  and  strain-rate  conditions: 
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If  p2  »  1.  then  Equation  (A2)  reduces  to 
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where 


For  the  limiting  condition  where  the  strain  and  therefore  M  are  large  (so 

that  tanh  M  =  I),  Equation  (A3)  becomes 
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where  Oj  ia  defined  as  the  limiting  atreaa  aaaociated  with  large  atraina. 
Finally,  if  Zp  »  1,  then 
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C  Creep 

For  the  condition  where  the  temperature  and  load  (stress)are 
constant,  Equation  (Al)  yields  the  following  relation  for  the  secondary, 
steady  creep  rate: 
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If  o-  »  (t0,  Equation  (A6)  reduces  to 

*«"*“- (f-i?) 
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or 
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(AS) 


Equation  (Al),  in  itself,  does  not  give  information  concerning  the 
final  or  tertiary  stage  of  creep  nor  the  creep  rupture.  However,  it  is 
possible,  with  the  help  of  some  qualitative  consideration,  to  arrive  at  an 
equation  for  creep  rupture.  For  example,  assuming  that  the  strain- rate 
in  the  final  stage  (4f)  is  of  the  form 
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where  p  is  a  constant  and  taking  the  function  to  be 


we  have 


(A10) 
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Integration  of  this  relation  yields 

*l ' exp  (' ?) 
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where  tf  is  the  tertiary  stage  time.  If  Cf  =  er  »  p  when  the  time  to  rupture 
is  tf  =  tr,  then  from  Equation  (A12),  we  have 

P  =  ectr  =  constant  (A13) 

Thus,  with  the  values  for  the  steady  creep  rate  determined  from  Equa¬ 
tions  (A6)  or  (A7)  and  the  extra  constant  p,  it  should  be  possible  to  calcu¬ 
late  the  creep  rupture  time. 

D.  Stress  Relaxation 

For  those  test  conditions  where  the  strain  and  temperature  are 
maintained  constant  and  the  stress  (load)  changes  as  a  function  of  time, 
we  have  from  Equation  (Al)  (for  e  =  0) 
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When  integrated,  this  yields  the  following  expression  for  the  rate  of  decay 


of  stress  with  time. 
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where  <r(t0)  is  the  applied  stress  at  t  =  0  and 
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E.  Increasing  Temperatures 

For  the  condition  where  the  load  (stress)  is  constant  and  the  temper¬ 
ature  increases  at  a  constant  rate  ('f  =  I'd),  Equation  (Al)  yields  the  following 


relation  for  the  strain  at  any  absolute  temperature  (T): 


(A17) 


where 


If  a  rod  is  restrained  at  the  end  and  is  subjected  to  a  linearly  increas¬ 
ing  temperature,  the  change  in  thermal  stress  (Acr)  may  be  acquired  in  the 
step-by-step  solution  of  the  following  expression: 

A,  =  E[a.(^)exp(X.|H)]iT  (A18) 

F.  Periodic  Stress 

For  the  condition  where  the  imposed  load  may  be  expressed  as  a 
periodic  function: 

<r  =  A  sin  wt 

When  used  in  Equation  (Al),  we  have 
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€  = 


u>  A 


cob  wt  +  2bT  exp  (  -■—[  )  ainh  (A  sin  wt)  (A19) 


Equation  (A19)  would  indicate  that  there  is  no  net  effect.  This  is  due  to  the 
equal  value  placed  on  the  positive  and  negative  parts  of  the  cycle.  During 
the  compression  portion  of  the  cycle,  however,  the  greater  part  of  the 
creep  during  the  tension  portion  is  not  recovered.  Consequently,  a  positive 
result  is  obtained,  if  the  contribution  of  the  compression  portion  is  consid¬ 
ered  as  zero.  The  strain  per  cycle  is  then 
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where  t  is  the  period  of  stress  alternation  and 


(A20) 


In  being  the  modified  Bessel  function. 

Assuming  rupture  occurs  in  the  same  manner  as  when  the  metal  is 
under  constant  stress: 

^clr  =  €r 

and  the  mean  value  of  ec  over  a  cycle  is 


and  also 

tr  =  Nrr 

where  Nr  is  the  number  of  cycles  required  for  rupture.  Therefore, 
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€c*r  =  cpc  =  er 
or 

cr=Nr  .Tt  e^,(.|H)F(A)  <A2J> 

This  equation  contains  only  known  quantities  with  the  exception  of  Nr  and, 
theoretically,  could  be  solved  for  Nr.  Thus,  for  the  first  time  we  have 
a  fatigue  equation  with  no  empirical  factors. 
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